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ABSTRACT 


The  Van  Atta  reflector  was  first  described  in  a  patent  by  Dr.  L.C. 

Van  Atta  in  195 9«  The  advantage  of  this  passive  reflector  type  should  be 
that  the  reradiated  field  has  a  maximum  back  in  the  direction  of  arrival 
of  the  primary  plane  wave.  Since  this  retrodirective  effect  of  the  reflec¬ 
tor  might  bo  of  great  importance  if  used  os  a  navigational  aiu  in  the  air 
or  at  sea,  it  t>  med  worth  while  to  carry  out  a  theoretical  investigation 
of  such  reflectors,  especially  since  only  experimental  investigations  had 
been  made  before  this  contract  was  initiated. 

The  work  performed  under  the  contract  deals  mainly  with  theoretical 
and  numerical  investigations  of  Van  Atta  reflectors  consisting  of  dipoles. 

A  survey  of  the  literature  concerning  active  or  passive  Van  Atta  reflectors 
has  been  made.  Both  a  linear  and  a  two-dimensional  plane  Van  Atta  reflec¬ 
tor  has  been  investigated  numerically  and  a  theory  for  arbitrary  Van  Atta 
reflectors  has  been  developed.  An  experimental  investigation  of  a  linear 
Van  Atta  reflector  was  carried  out  and  the  results  compared  with  the  theo¬ 
retical  results. 
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1.  INTRODUCTION 

The  contract  AP  6l(052)-79*»  wan  running  from  April  ll(  196b  to  June 
30,  196T. 

The  objectives  o+  the  contract  were  given  as  follows: 

Carry  out  a  theoretical  investigation  of  the  performance  of  a  Van 
Atta  reflector  and, by  electronic  computer  calculations,  to  find  the  radia¬ 
tion  patterns  of  a  number  of  such  reflectors  with  various  characteristics. 
It  is  expected  that  this  would  lead  to  a  procedure  for  the  design  of  a  Van 
Atta  reflector  with  a  prescribed  radiation  pattern.  It  is  the  intention 
first  to  treat  the  linear  and  two-dimensional  array,  later  the  method  of 
investigation  may  be  extended  to  include  Van  Atta  refleotors  placed  on  cir¬ 
cular  cylinders  and  spheres.  Further  efforts  will  be  i..ade  to  supplement 
the  theoretical  investigation  by  experiment. 

In  this  final  report  a  summary  of  all  the  work  performed  under  the 
contract  is  given.  Some  related  work  carried  out  at  this  laboratory  but 
not  paid  by  the  contract  is  described,  too. 

First  an  analytical  investigation  of  arbitrary  Van  Atta  reflectors 
was  carried  out.  Each  pair  of  antenna  elements  with  connecting  transmis¬ 
sion  lines  was  described  by  an  equivalent  circuit.  Reradiation  from  the 
elements  and  mutual  interaction  were  taken  into  account.  Next  a  theoreti¬ 
cal  investigation  of  a  linear  Van  Atta  reflector  consisting  of  four  paral¬ 
lel  half-wave  dipoles  was  performed.  The  dipole  elements  were  equispaced 
and  mutual  impedances  between  the  dipoles  were  taken  into  account.  An 
expression  for  the  reflected  field  was  set  up  by  superposition  of  the  three 
fields  involved.  Further  an  analysis  of  the  general  shape  of  the  reradia?- 
tion  pattern  was  xaade  neglecting  the  mutual  impedances.  This  analysis 
showed  that  for  some  combinations  of  the  parameters  involved  the  reflector 
does  not  act  as  stated  in  the  patent  description  of  Van  Atta  reflectors. 

An  experimental  investigation  of  the  four  element  linear  Van  Atta 
reflector  was  carried  out  in  a  radio-anechoic  chamber.  In  agreement  with 
the  theoretical  investigation,  the  experiments  showed  that  the  reflector 
only  to  a  limited  extent  has  the  retrodirective  effect  stated  in. the  patent 
description  and  that  it  has  a  mirror  effect  to  the  sane  extent  as  it  has  r. 
retrodirective  effect. 

Numerical  investigations  of  a  four  element  Van  Atta  reflector  and  an 
optimisation  of  the  reflector  with  respect  to  the  parameters  involved  was 
carried  out  using  an  electronic  computer.  The  reradiation  pattern  of  the 
reflector  vaB  optimized  with  the  criterium  that  the  minimum  values  of  the 
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back-scattered  field  intensities  for  all  angles  of  incidence  should  be  as 
large  as  possible. 

Finally  a  square  Van  Atta  reflector  consisting  of  half-wave  dipoles 
was  investigated  theoretically  and  numerically.  The  effect  of  mounting 
the  dipoles  above  and  parallel  to  a  conducting  plate  was  also  examined. 

The  mutual  impedance  between  the  dipole  elements  and  the  reradiation  both 
from  the  elements  and  from  the  plate  was  taken  into  account.  All  parame¬ 
ters  involved  have  been  varied  in  a  numerical  analysis  of  the  reflector. 

The  parameters  are;  the  number  of  dipole  elements,  the  length  of  the 
dipoles,  the  length  and  characteristic  impedance  of  the  transmission  lines, 
the  distance  between  the  dipoles,  and  the  distance  from  the  dipoles  to  the 
conducting  plate.  The  effect  of  changing  the  parameters  is  shown  in  curves 
of  the  back-scattering  cross  section  of  the  reflector  as  a  function  of  the 
variation  of  each  of  the  parameters.  The  numerical  results  have  been 
compared  with  experimental  results  obtained  by  others. 

The  last  work  cn  the  contract  has  been  an  investigation  of  the  band¬ 
width  properties  of  a  k  by  k  element  square  Van  Atta  reflector  consisting 
of  dipoles. 


2.  SURVEY  OF  LITERATURE 


8inco  Van  Atta  ^  in  1959  proposed  bis  passive  retrodirective  reflector 
several  papers  have  suggested  the  use  of  this  reflector  type  in  various 
caamunication  systems  2  23)^  8cmo  of  the  papers  tend  to  give  an  analy¬ 

tical  treatment  of  the  Van  Atta  reflector  but  most  of  them  neglect  the 
scattering  by  and  coupling  between  the  antennas «  Many  of  the  papers 
suggest  active  component a  inserted  in  the  transmission  lines  of  the  reflec¬ 
tor,  such  as  modulated  phase  shifters  5\  amplifiers  ^  10 \  and  mechanical 
modulation  by  means  of  cavity  resonators  It  has  bean  suggested  to  use 
Van  Atta  reflectors  for  satellite  communication  and  both  passive,  semlpas- 
sive  and  active  systems  have  been  proposed  V  15 *  21  *  2Z\  Other  applica¬ 
tions  are  for  navigational  aids,  for  example  used  to  enhance  the  reflec¬ 
tion  from  radartargets,  from  lifeboats,  and  from  aircraft  10^.  A 

more  detailed  discussion,  of  the  literature  concerning  Van  Atta  reflectors 
and  their  applications  is  given  in  Scientific  Report  Ho.  1  and  by 
Appel-Eansen  ^*0  ^ , 
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3.  THEORETICAL  INVESTIGATION  OF 
ARBITRARY  VAN  ATTA  ARRAYS 


A  general  treatment  which  nay  be  uaed  for  a  number  of  different 
investigations  is  described  in  Scientific  Report  No.  1  of  the  contract. 

The  configuration  investigated  is  shown  in  fig.  1.  The  elements  are  sup¬ 
posed  to  be  dipoles ,  but  the  theory  could  easily  be  extended  to  other  an¬ 
tenna  types,  the  dipoles  arc  placed  on  and  parallel  to  an  imaginary  smooth 
surface  which  may  be  f.ex.  a  plane,  a  cylinder  or  a  sphere.  The  field 
incident  on  the  reflector  is  a  plane  wave. 

The  open  circuit  voltage  induced  at  the  terminals  of  each  antenna 
element  by  the  primary  plane  wave  is  calculated.  Using  these  voltages  a 
system  of  linear  equations  is  developed  for  calculating  the  currents  in  each 
antenna  taking  into  account  the  mutual  impedances,  the  characteristics  of 
the  interconnecting  transmission  lines,  and  the  induced  voltage  at  the  element 
itself  (giving  the  scattered  field)  and  at  its  mate  (giving  the  retrodirective 
field). 

When  the  currents  are  determined  the  reradiated  field  may  be  calculated. 
For  a  reflector  vith  all  elements  parallel  this  is  done  using  the  theory  of 
antenna  arrays.  Finally  the  properties  of  the  reflector  array  are  described 
by  calculation  of  the  differential  scattering  cross  section. 

The  mutual  impedances  of  the  dipoles  are  calculated  using  the  induced 
EMF  method  as  given  in  Jordan's  book  .  Algol  procedures  for  computing 
the  mutual  impedances  between  linear  dipoles  vith  sinusoidal  current  distri¬ 
butions  and  for  arbitrary  wire-antennas  vith  a  known  current  diatribution 
has  been  developed  at  this  laboratory. 

The  transmission  lines  are  represented  by  equivalent  circuits  of  the 
general  x-circuit  type.  Ibis  type  of  equivalent  circuit  has  been  chosen 
because  it  has  the  advantage  of  being  valid  for  all  lengths  of  the  trans¬ 
mission  lines.  It  is  assumed  to  be  symmetrical  and  lossless. 

Since  both  the  induced  voltages  and  the  mutual  impedances  are  complex 
quantities  the  real  matrix  equation  to  be  solved  will  be  of  the  2Rth  order 
where  N  is  the  total  number  of  antenna  elements.  Ibis  means  that  it  vill  be 
almost  necessary  to  use  an  electronic  computer  for  solving  the  matrix  equa¬ 
tion  numerically  vhen  reflectors  vith  more  than  tvo  elements  are  treated. 

A  simple  numerical  example  of  a  four  element  Van  Atta  array  vith  equi- 
spaed  half-wave  dipoles  is  given  in  SR  1  in  order  to  illustrate  the  theory 
developed.  The  numerical  results  indicate  that  a  retrodirective  effect  as 
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stated  by  Van  Atta  is  obtained  to  acme  degree.  However  the  results  also 
show  that  the  mirror  effect  of  the  reflector  ib  of  the  some  order  of  magni¬ 
tude  am  the  retrodirectivc  effect.  The  influence  of  the  mutual  impedances 
and  of  a  mismatch  between  the  antenna  elements  and  the  transmission  lines 
may  be  utilized  to  change  the  reradiation  pattern  to  compare  better  with  a 
prescribed  form. 


h.  A  LINEAR  VAN  ATTA  REFLECTOR 
CONSISTING  OF  FOUR  HALF-WAVE  DIPOLES 


A  detailed  investication  of  a  four  element  linear  Van  Atta  reflector 
consisting  of  half-vave  dipoles  was  carried  out  and  is  deacribed  in  Scien¬ 
tific  Report  IJos.  2,  3,  and  L  of  the  contract. 

Fij  st  an  expression  for  the  reflected  field  is  devoloped  considering 
each  dipole  as  a  receiving  and  transmitting  antenna  matched  to  its  trans¬ 
mission  line.  When  a  plane  wave  is  incident  on  the  reflector  a  current 
distribution  consisting  of  three  different  parts  will  be  generated  in  the 
antenna  elements.  The  first  part  is  due  to  the  energy  transmitted  to  the 
antenna  through  the  transmission  line  from  itB  mate.  The  second  part  is 
the  current  induced  in  the  antenna  by  the  incident  plane  wave  and  the  third 
part  is  induced  in  the  antenna  due  to  its  mutual  interaction  with  the  other 
antennas  of  the  reflector.  Only  the  first  part  of  the  current  distribution 
is  considered  in  Van  Atta's  patent  description,  while  the  last  two  parts  are 
neglected. 

In  fig.  2  is  shown  the  retrodirective  effect  of  the  reflector  due  to 
the  first  part  of  the  current  distribution  mentioned  above.  The  second  part 
of  the  current  distribution  creates  a  mirror  effect  of  the  reflector  which 
when  the  dipoles  are  matched  to  the  transmission  lines  is  of  the  same  order 
of  magnitude  as  the  retrodirective  effect.  This  mirror  effect  which  is 
shown  in  fig.  3  is  not  mentioned  in  the  patent  description  of  the  reflector. 

It  is  shown  that  the  length  of  the  transmission  lines  and  the  angle 
of  incidence  has  a  great  influence  on  the  rcradiation,  and  for  some  values 
of  these  tvo  parameters  the  first  and  second  part  of  the  current  in  the 
antennas  are  of  opposite  phase  and  they  cancel  each  other  so  that  the  only 
reflected  energy  is  the  small  part  due  to  the  mutual  interaction  between  the 
elements.  The  length  of  the  transmission  lines  at  which  the  behaviour  of 
the  reflector  is  as  much  as  possible  in  accordance  with  the  patent  description 
has  been  found  but  even  this  reflector  has  the  mirror  effect  mentioned  above. 

Numerical  analysis  have  been  made  with  □  spacing  between  the  dipole 
elements  of  half  a  wavelength.  It  turns  out  that  in  most  cases  maximum 
reflection  is  not  back  in  the  direction  of  incidence  of  the  primary  plane 
wave,  and  that  the  mutual  interaction  between  the  dipoles  causes  asymme¬ 
tries  in  the  reradiation  pattern.  This  means  that  it  is  possible  to  increase 
the  back-scattered  energy  by  choosing  a  proper  combination  of  tho  distance 
between  the  elements  and  the  length  of  the  transmission  lines.  The  above 


mentioned  investigations  are  described  in  Scientific  Report  Ho.  2. 

In  order  to  verify  the  theoretical  results  obtained, an  experimental 
investigation  of  a  four  element  linear  Van  Atta  reflector  was  carried  out 
in  a  radioanechoic  box  at  the  laboratory.  This  investigation  is  treated 
in  Scientific  Report  No.  3.  The  four  half-wave  dipoles  of  the  experimental 
reflector  were  slot  fed  dipoles  with  open-ended  terminations.  The  length 
of  the  connecting  transmission  lines  could  be  changed  by  means  of  line- 
stretchers  thus  examining  the  influence  of  the  length  of  the  lines  on  the 
reflecting  properties  of  the  reflector.  The  measurements  were  performed  at 
3.21  GHz  as  this  frequency  gives  the  best  matching  between  each  dipole  and 
the  connected  transmission  line. 

The  Van  Atta  reflector  was  placed  on  a  moveable  pedestal  in  the  center 
of  the  anechoic  box,  and  the  measurements  were  based  upon  the  principle  of 
interference  between  the  signal  reflected  by  the  reflector  and  a  reference 
signal.  Radiation  patterns  were  measured  in  the  plane  normal  to  the  axis 
of  the  dipoles  for  descrete  values  of  the  angle  of  incidence  of  the  primary 
plane  vave. 

A  good  agreement  between  the  experimental  and  theoretical  results  was 
found.  The  measured  results  confirmed  the  theoretical  results  of  Scientific 
Report  Ho.  2  and  showed  (l)  that  maximum  reradiation  is  not  always  back  in 
the  direction  of  incidence,  (2)  that  the  reflector  has  a  mirror  effect  to 
the  same  degree  as  it  has  a  retrodirective  effect,  (3)  that  the  reradiation 
depends  strongly  upon  the  length  of  the  transmission  lines,  and  (U)  that  the 
mutual  impedances  causes  asymmetries  in  'he  radiation  patterns. 

As  an  extension  of  the  theoretical  analysis  of  the  four  element  linear 
Von  Atta  reflector  an  optimization  of  this  reflector  has  been  carried  out, 
the  results  of  which  is  given  in  Scientific  Report  Ho. 

First  the  original  expression  for  the  reradiated  field  as  derived  in 
Scientific  Report  No.  2  was  changed  to  a  more  general  form  which  makes  it. 
possible  to  study  the  influence  of  asymmetries  in  the  location  of  the  Clipo  ves , 
unequal  lengths  of  the  transmission  lines,  and  a  mismatch  between  the  dipc.ef 
and  the  transmission  lines.  Further,  a  method  was  developed  for  computing 
a  quantity  which  may  be  used  as  a  measure  of  the  deviation  from  the  retro- 
directive  effect  of  the  reflector.  This  quantity  is  shown  to  be  just  as 
useful  as  the  reradiation  pattern  itself  when  two  different  Van  Atta  reflec¬ 
tors  are  to  be  compared, 

A  perfectly  working  retrodirective  Van  Atta  reflector  has  not  been 
found.  However,  when  mutual  coupling  is  neglected,  a  condition  for  the 
smallest  deviation  from  retrodirectivity  has  been  derived.  It  i3  shown  that 
the  mutual  coupling  between  the  dipoles  usually  causes  the  reradiation  of 
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the  reflector  to  decrease  and  the  deviation  from  retrodirectivity  to  increase. 
However  for  certain  values  of  the  parameters  involved  it  turns  out  that  coupl¬ 
ing  may  increase  the  back-scattcring  up  to  50,1  for  some  ancles  of  incidence. 

The  numerical  optimization  of  the  reradiation  pattern  of  the  four 
element  Van  Atta  reflector  was  carried  out  using  a  computational  technique 
starting  with  on  a  priori  reasonable  set  of  parameters  selected  by  examininc 
1600  different  reflectors.  The  parameters  were  then  perturbed  about  their 
initial  values,  the  effect  on  the  reradiation  pattern  was  observed,  and  a  new 
set  of  parameters  giving  an  improved  result  was  selected.  The  success  of  this 
method  depends  cm  the  correctness  of  the  original  set  of  parameters  and  the 
computer  program  for  perturbing  the  parameter  values 

An  attempt  was  made  to  fulfil  the  following  two  criteria: 

the  minimum  value  of  the  back-scattered  field  intensity,  as  a  function 

of  the  angle  of  incidence,  should  be  os  large  as  possible  and  the 

deviation  from  Van  Atta  effect  as  small  as  possible, 

the  minimum  value  of  the  back-scattered  field  intensity,  as  a  function 

of  the  angle  of  incidence,  should  be  above  various  prescribed  levels 

and  the  deviation  from  Von  Atta  effect  as  small  as  possible. 

For  both  optimization  processes  it  turned  out  that  the  optimum  value 
of  the  spacing  was  close  to  1.5  wavelengths.  Further  it  was  found  that,  due 
to  coupling,  the  minimum  value  of  the  back-scattered  field  intensity  may  be 
increased  and  the  deviation  from  the  retrodirective  effect  decreased  if  the 
transmission  lines  are  permitted  to  be  of  unequal  lengths  and  asymmetries 
are  permitted  in  the  location  of  the  dipoles  around  the  center  of  the  reflec¬ 
tor.  However,  the  improvements  are  snail  and  asymmetries  often  causes  the 
opposite  effect. 
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5.  SQUARE  VAN  ATTA  RKFLECTORG  I 

WITH  OR  WITHOUT  A  COIiDUCTIHG  PLATE 

1 

Another  configuration  of  the  reflector  array  is  the  plane,  square 
Van  Atta  reflector  consisting  of  parallel  dipoles.  The  investigation  of 
this  reflector  is  described  in  Scientific  Report  No.  5. 

The  theoretical  investigation  of  this  reflector  has  already  been  j 

described  in  Scientific  Report  No.  1.  However,  in  order  to  be  able  to  j 

2)  ! 

compare  the  theory  with  experimental  results  obtained  by  Sharp  the  ! 

effect  of  mounting  the  dipoles  above  and  parallel  to  a  conducting  plate 
has  to  be  taken  into  account.  The  system  investigated  is  shown  in  fig.  I*. 

the  reflecting  properties  of  the  plate  are  supposed  not  to  be  in¬ 
fluenced  by  the  presence  of  the  dipoles.  The  reflected  field  is  found 

•  39 ) 

using  the  method  of  physical  optics  as  described  f.ex.  in  Kerr's  book 

for  a  plate  the  dimensions  of  which  are  not  small  compared  to  the  wave¬ 
length. 

The  reradiating  properties  of  the  dipoles  when  the  plate  is  present, 
is  calculated  as  if  the  plate  was  infinite  in  extent,  using  the  theory  of 
images.  The  system  of  dipoles  may  then  be  treated  along  the  same  lines  a>: 
in  Scientific  Report  No.  1,  but  the  induced  voltage,  the  mutual  impedance:  , 
and  the  determination  of  the  field  reradiated  from  the  dipoles  have  to  be 
changed  because  of  the  image. 

The  induced  voltage  is  still  found  as  described  for  an  arbitrary  j 

reflector  in  Scientific  Report  Ho.  1,  but  now  the  electromagnetic  field  , 

vector  is  changed  in  such  a  way  that  the  distance  from  the  dipoles  to  the 
plate  is  involved,  according  to  ordinary  reflection  theory.  j 

The  new  values  of  self-  and  mutual  impedances  are  found  using  the  method  ! 

of  images,  too.  *  j 

By  using  the  values  of  the  induced  voltages  and  the  values  of  the  s:  lf- 
and  mutual  impedances  thus  found,  the  system  of  equations  (23)  of  Scientific  , 

Report  No.  1  will  give  the  new  values  of  the  currents  on  the  arttenna  elements  j 

when  the  presence  of  the  plate  is  taken  into  account.  ! 

After  that  the  reradiation  pattern  from  the  dipole  reflector  itself 
mounted  in  a  distance  li  above  the  plate  is  calculated  vith  the  above-mentioned 

currents  on  the  dipole6.  The  final  reradiation  pattern  of  the  dipoles  is 

,  i 

found  using  the  theory  of  antenna  arrays  on  the  array  consisting  of  two  paral¬ 
lel  Van  Atta  reflectors  in  free  space  with  the  distance  2h,  where  h  is  the  j 
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distance  between  the  dipoles  and  the  plate. 

The  total  field  reradiated  from  the  reflector  system  is  found  by  , 

adding  the  field  reradiated  from  the  dipoles  and  the  field  reflected  from 
the  conducting  plate. 

Using  the  above-mentioned  theory  a  computer  program  has  been  developed 
and  the  numerical  results  have  been  compared  vrith  results  obtained  by  Sharp 
from  experimental  investigations  of  a  16  element  square  Van  Atta  reflector. 

The  computed  back-scattering  cross  section  shows  a  good  agreement  with  the 
results  measured  for  the  experimental  reflector  as  shown  in  fic.  5. 

Furthermore,  a  series  Of  computations  has  been  performed  in  order  to 

examine  the  changes  in  the  back-scattering  cross  section  due  to  changing  of  j 

j 

the  parameters  of  the  reflector.  The  parameters  are  the  number  (N)  of  j 

elements,  the  length  (a)  and  characteristic  impedance  (Zq)  of  the  trans-  j 

mission  lines,  the  distance  (d)  between  adjacent  dipole  elements,  and  the 
distance  (h)  from  the  elements  to  the  plate. 

The  most  important  results  obtained  by  this  numerical  investigation  is: 

(1)  that  the  back-scattering  cross  section  becomes  larger  if  a 
distance  of  0.35  wavelengths  from  the  dipoles  to  the  plate  is 

used  instead  of  0.25  wavelengths  as  used  by  Sharp,  j 

(2)  that  the  shape  of  the  curves  of  back-scattering  becomes  more  j 

| 

irregular  when  more  elements  are  used  in  the  reflector  but  the 
level  of  back-scattering  i3  increased, 

(3)  that  the  shape  of  the  curves  of  back-scatterinc  becomes  more 

smooth  when  a  mismatch  between  the  elements  and  the  transmission  j 

lines  is  introduced  in  such  a  way  that  the  characteristic  impe¬ 
dance  of  the  line  is  larger  than  the  self~impedar.ee  of  the  dipole. 

However,  the  magnitude  of  the  back-scattered  energy  is  then  de¬ 
creasing  more  rapidly  for  oblique  directions  of  incidence* 

(h)  that  for  certain  lengths  of  the  transmission  lines  the  back- 

scattering  in  the  direction  normal  to  the  reflector  tends  to  zero. 

This  is  in  accordance  with  the  results  obtained  for  the  four 
element  linear  reflector  mentioned  in  section  1». 

In  Scientific  Peport  Ho.  5  a  great  number  of  numerical  results  obtained 
by  the  parameter  variation  is  given  as  curves  of  the  back-scattering  cross 
section. 
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6,  DAHWIDTH  PROPERTIES  OF  THE  SQUARE 
VAH  ATTA  REFLECTOR 


The  final  work  on  the  contract  deals  with  a  computation  of  the  band- 
vidth  properties  of  a  l6  element  square  Van  Atta  reflector  similar  to  the 
experimental  model  used  by  Sharp. 

Part  of  this  work  is  on  investigation  of  the  influence  on  the  back- 
scattcrinc  properties  of  the  reflector  of  a  change  in  the  length  of  the 
dipole  elements.  In  the  previous  investigations  on  this  contract  only 
reflectors  consisting  of  half-wave  dipoles  have  been  investigated.  The 
examination  of  the  influence  of  the  length  of  dipoles  is  carried  out  along 
the  same  lines  as  the  examination  of  the  influence  of  the  other  parameters 
of  the  reflector  as  described  in  Scientific  Report  Mo.  5.  Thi3  means  that 
the  length  of  the  dipoles  is  changed  while  keeping  all  other  parameters  of 
the  reflector  fixed  with  values  corresponding  to  the  dimensions  of  the 
experimental  reflector  used  by  Sharp. 

The  results  of  this  investigation  is  shown  in  fig.  6  of  this  report. 
From  this  it  turns  out  as  expected  that  the  back-scattered  energy  decrease ; 
both  when  the  dipole  length  is  less  than  and  greater  than  half  a  wavelength. 
This  is  due  to  the  fact  that  the  matched  half-wave  dipole  ha3  optimal  reru- 
diating  properties.  However,  when  the  dipoles  are  less  than  half  a  wave¬ 
length  the  mutual  coupling  between  the  dipoles  decrease  and  a  better  retro- 
directive  effect  of  the  reflector  is  obtained.  This  effect  is  further 
strengthened  because  of  the  mismatch  between  the  dipole  and  the  transmission 
line,  the  self  resistance  PA  of  the  dipole  being  less  than  the  characteristic 

innedance  Z  of  the  line.  This  is  in  accordance  with  the  results  of  the 
o 

investigation  of  the  four  element  linear  reflector  where  it  was  found  that 
the  deviation  from  retrodirective  effect  decreases  when  the  factor  R^/Z^ 
decreases. 

The  bandwidth  properties  of  the  16  element  square  Von  Atta  reflector 
similar  to  the  experimental  reflector  used  by  Sharp  with  or  without  the 
conducting  plate  has  been  computed  and  the  results  are  shown  in  figs.  7a 
and  8a.  The  curves  show  the  back-scattered  energy  for  different  angles 
of  incidence  as  a  function  of  1A0»  where  is  the  wavelength  corresponding 
to  the  center  frequency. 

For  increasing  values  of  A/A^  larger  than  1.0  the  back-scattered  energy 
from  the  reflector  without  a  conducting  plate  decreases  in  a  regular  manner 
for  all  angles  of  incidence. 


For  values  of  X/X^  <  l.o  the  shape  of  the  curves  is  very  irregular 
and  give6  no  information  at  all.  When  the  conducting  plate  is  taken  into 
account  the  curves  for  A/Xq  >  1.0  ore  almost  as  regular  as  when  the  plate 
is  not  present  but  the  level  of  back-scattered  energy  is  higher  for 
directions  of  incidence  near  normal  incidence  (0°)  and  lower  for  directions 
of  incidence  near  broadside  (90°).  For  X/X^  <  1.0  the  curves  are  just  as 
irregular  as  in  the  case  where  the  plate  is  not  present. 

However,  from  the  results  of  this  investigation  it  turned  out  that 
the  retrodirective  effect  cf  this  16  clement  reflector  is  essentially  im¬ 
proved  when  X/Xo  «  1.3.  This  corresponds  to  a  16  element  square  Van  Atta 
reflector  with  the  following  parameter  values: 


length  of  dipoles 

0.3C5 

X 

radius  of  dipoles 

0.0115 

X 

distance  between  dipoles 

0.fc62 

X 

distance  from  dipoles  to  plate 

0.192 

X 

length  of  transmission  lines 

0.315 

X 

characteristic  impedance  of 

transmission  lines 

73.0 

ohms 

A  measure  of  retrodirective  effect  is  given  in  figs.  7b  and  Ob  as 

'l  =  number  of  examined  angles  of  incidence  riving  retrodirective  effect  § 

total  number  of  angles  ot  incidence  examined  ”  *  X  " 

For  the  reflector  without  a  plate  II  obtaines  its  maximum  90%  for  the  above- 
mentioned  wavelength  X  =  1. 3Xq  of  the  incident  plane  wave.  When  the  con¬ 
ducting  plate  is  present  the  optimum  value  of  M  is  50%  and  this  value  is 

obtained  in  the  whole  range  fro®  X  -  1.1X  to  X  «=  1.5X  .  The  reason  for 

o  o 

the  decrement  of  M  in  the  second  case  is  that  the  retrodirective  effect  is 
reduced  when  the  direction  of  incidence  turns  towards  broadside  because  of 
interference  with  the  field  scattered  from  the  conducting  plate. 

However  it  is  obvious  in  both  cases  that  a  Van  Atta  reflector  with  a 
better  retrodirective  effect  than  the  effect  measured  by  Sharp  may  be 
obtained  from  the  same  physical  reflector  at  the  expense  of  the  reradiated 
energy  if  the  reflector  is  used  at  lower  frequencies  than  the  center-fre¬ 
quency. 
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7.  CONCLUSION 

In  this  report  a  survey  of  the  investigations  performed  under 
Contract  AF  61(052) -T91*  "Reflector  Array"  has  been  given  including  the 
results  of  the  final  vork  on  the  contract  vhioh  have  not  been  described 
in  any  previous  report. 

By  comparing  the  results  obtained  with  the  objectives  of  the  con¬ 
tract  as  stated  in  section  1  of  this  report  it  is  seen  that  the  first 
period  of  the  text  may  be  covered  by  Scientific  Reports  Nos.  1.  2.  4,  and 
5.  The  second  period  is  in  part  covered  by  Scientific  Report  No.  4.  The 
linear  and  two-dimensional  array  mentioned  in  the  third  period  has  been 
investigated  while  the  other  configurations  mentioned  in  this  period  have 
not  been  dealt  with.  The  fourth  and  last  period  is  covered  in  Scientific 
Report  -No.  3  for  a  linear  array. 

Further  an  investigation  of  the  bandwidth  properties  of  a  square 
reflector  has  been  carried  out. 

Using  the  theory  explained  in  the  reports  issued  under  this  contract 
it  should  be  possible  to  investigate  other  types  of  Van  Atta  reflectors. 
Ibis  might  be  as  well  reflectors  consisting  of  dipoles  mounted  in  two-  or 
three-dimensional  arrays  for  example  over  conducting  cylinders  or  spheres, 
as  reflectors  consisting  of  ocher  types  of  antenna  elements  such  as  horns, 
paraboloid  antennas,  crossed  dipoles  or  monopoles. 

Probably  results  which  compare  better,  with  the  experimental  results 
measured  by  Sharp  may  he  obtained  using  another  theory  for  the  field  re¬ 
flected  from  the  conducting  plate  than  the  physical  optics  theory  used  in 
Scientific  Report  No.  5«  This  theory  may  be  Keller's  geometrical  theory  of 
optics  which,  in  contrast  to  the  theory  used,  will  take  into  account  the 
scattering  of  the  incident  field  about  the  edges  of  the  conducting  plate. 
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8.  COMPUTER  PROGRAM 


In  Appendix  1  &  copy  of  the  computer  pro eras  developed  for  the  nume¬ 
rical  investigation  of  square  Van  Atta  reflectors  vith  or  vithout  a  con¬ 
ducting  plate  is  printed. 

The  program  is  written  in  FORTRAN  IV  and  the  computer  used  in  an 
IBM  7090  run  fcy  the  Northern  Europe  University  Computing  Center,  Techni¬ 
cal  University  of  Denmark. 
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COB  *  STEPS  IN  CHANGING  8  DURING  PARAMETER  VARlATICN 

02  =  STEPS  IN  CHANGING  ZC  OURING  PARAMETER  VARlATICN 

CCLE  =  STEPS  IN  CHANGING  DIE  DURING  PARAMETER  VAR  I  AT  I Of 

IMX  =  MAXIMUM  VALUE  CP  I  IN  PARAMETER  VARIATION 

CMX  =  MAXIMUM  VALUE  OF  D  IN  PARAMETER  VARIATION 

RMX  =  MAXIMUM  VALUE  OF  R  IN  PARAMETER  VARlATICN 

BMX  =  MAXIMUM  VALUE  OF  B  IN  PARAMETER  VARlATICN 

ZMX  =  MAXIMUM  VALUE  OF  ZC  IN  PARAMETER  VARlATICN 

CLMX  =  MAXIMUM  VALUE  OF  DUE  IN  PARAMETER  VARIATION 


DIMENSION  A( 72,72) ,C( 72) ,X (72) ,U (72,10) 

DIMENSION  H( 3f ,6) ,Q(2, 1C>  »S(5) ,V(36,10) ,ie (6) 

DIMENSION  G(  IS,  10)  t  TETAI  ( 1C )  .  PH  1 1  ( 10 )  ,  VV  ( 10  )  ,T  1 19 )  ,  F  (  »,  9  \ 
INTEGER  BW , C I , PAG  I NA 
369  FORMAT (315) 

301  FORMAT  <  A 1 5  > 

302  FORMAT ( 4F5 • 3 ) 

303  FORMAT ( 2F9  • 3 ) 

304  FORMAT ( 3F9 ♦ 3  ) 

305  FORMAT! 6FS. 3) 

306  FORMAT (15,  5F9, 3  I 

ICO  REAC(5,3C1)  I , J»KPQ,MDB 

READ! 5, 3C5  )  D,R»B, ZO,RAC  ,DLE 
REAC( 5» 3C3  )  W,FII 
READ ( 5»  369 )  LFI »MA,BW 
C 

IF(BW.EQ.l)  READ ( 5 , 301 )  MDI P , MD , NR  ,MB 
IFIdW.EQ.l!  READ! 5, 3C3 )  DIFF, RATIO 

0 

IF( 3W.EQ.C )  READ ( 5  , 3C6 )  DI ,CD »DR ,CCB ,DZ  ,CC L E 
IF( BW.EQ  .C  I  READ! 5»3Cfc)  IMX, OMX, RNX»BNX , ZMX, CLMX 

IF(BW.EQ.C)  DIFF  =  0. 

IF(BW.EQ.C)  MC IP  =  1 
K=  IC 

LTETA*10 

FAKT0R=1.C-DIFF 

PAGINA  =  C 

I  ST  =  I 

CST  =  C 

RST  =  R 

BST  =  B 

ZST^ZO 

CLST  =  CI.  E 

c 

IF(DW.F0-1J  CO  TO  ‘3  02 
IF(BW.EQ.C)  GO  7C  5C3 
501  1=1  +CI 
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GO  TO  515 

510  C=C  +  CD 
GO  TO  515 

512  e=e  +cob 

GO  TO  515 

513  ZO=Z04DZ 
GO  TO  515 

511  R=R  4PR 
GO  TC  515 

51 A  CLE  =  CLE4CCLF 

c. 

515  IF(DW.EQ.C)  GC  TO  5C3 

502  IF(MC.EQ.C)  C*C/FAKTCR 
IF(WR.EO.C)  RsR/FAKTCR 
IF(MB.EQ.C)  B*B/FAKTQR 

503  C=C«£. 2831853 
R=R«fc. 2831853 
B=B«6. 2831853 
PAGINA  =  PAGINA  +  1 
N= I «  *  2 

M=I/2 

IF(M*2-I)  11,12,11 

11  LE=1 

GO  TO  2 

12  LE  =  C 

2  MM  =  (  I-LE  ) /2+5/ ICC 
CO  13  M=1»FM 

13  $ ( M  ) »(  FLOAT  m/2.-FLCAT(M)  +  C.S)*n 
LL*( l-LE+ll/2 

CO  15  L=1,LL 
CO  15  M* 1 , 1 
L 1  =  1 L“ 1 )  *  I  +M 
ML  1,2)=  S  f  L  I 
12= I-L ♦ 1* ( M- 1 ) *  I 
HL2,  1  )=  S ( L  ) 

L3=N-I+M-tL-l)»I 
F(  L  3 »  2  )  =-$(  L  ) 

LA=L4(M-1)#I 
MLA,  1 )  =- S  (  L  ) 

IF(LE)  1  A, 15, 1A 
1A  L5=(  I-1)/24(M-1)*I  +  1 
HL5,  1  )=C.C 
L6=(  IM 1-11 1/24M 
M L£, 2)=C.C 
15  CONTINUE 
I  I=I»«2 

FAKT1  =  l.CCCCCC 

IF(MCIP.EQ.O)  FAKT1  *  FAKTCR 

CALL  EETAUI, J,ZC,R,H,FAKT1  ,RAO,CLE) 

TAL=0. 

0PG=FL0AT(KPQ)/1C. 

TAL=TAL* 1. 

CLE=CLE/FAKT1 
CO  2 A  W=1,K 

YM  =  y 
VV( H >*W 


IFMS) 


no  n 
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PHIIIM»«FII 
XM  »  YM-dPQ 
XK=K- 1 

TETAi(MI*<XH*5C.)/XK 
IFUBStTETAKMIJ-SC.)  22,21,22 

21  IF(A8S(PHII(M>)— 9C.)  22,25,22 

22  TETIaTETAIt  M  )• C. Cl 7A  532925 
Z=0.C174532925 

OP=(COS(3«14159265*DLE*SINCTE7I)«SIN(FII«ZI)“CCS  { 3»14159265«CL  E  ) )  * 

1  (  SIN(W*Z  J«CCS(FII*Z)-CCS<W*2)*SIMFII#Z)«CCSnETl  M/ISCRTll.- 

2  SIN( TETI  1**2*SIN(F[ 1*Z)«*2)»S1K(3.14159265*CIE) ) 


Q(  1,M»=P«(  1. -FLOAT!  J)*CCS(2.«R*CCS(TCTID) 

Q( 2,M )a-p*FLCAT{ J)*SIN(2.*R*CCS (7E7I) ) 

CO  23  L*1 » N 

23  CV(L,M)s~F(L,l)»COS(FII«ZJ*SIN{TE7I) 

1  -H(1,2)«SIN(FII«Z)»S1MTETI) 

2A  CONTINUE 
GO  TO  26 

25  K*K- 1 

now  the  auxiliary  array  has  been  calculated 

26  CALL  MATKlEtltN, J,B,IB,A,H) 

IF(LG-l)  11C, 41, 11C 

41  CALL  MAT2(LE,I,N,J,B,A,H> 

110  CONTINUE 

C 

DO  50  L=  1, K 
1FILE-1 )  49,4E,49 
48  Nl*(N+l)/2 

N2*(2»N+l)/2 

CIN1I  a(Q(  1,U*CCS(V(N1,L)  )-G(2fU«SIMV(M,L) )  )«CCSCB) 
C ( N  2  >  a(Q(l,L)»SIN(V(Nl,l) ) +  G ( 2 , L ) *CCS ( V (N1 , L ) > )*CCS(B) 
*9  MM=(N-LEJ/2 
Nl»2«N+L 
CO  8CC  H a  1 , NH 
NMl=N  +  l-M 
MN=N*M 

Pl=M»(N+LE)/2 
M2=M  +  ( 3*N  +  LE  1/2 

RVN  =  Q(1»L)«CCS(V(M,L))-C(2,L)*SIN{V(P,L)) 

CVM  =  0I1,L)*SIN(V!M,L))*C(2,U*CCS(V(P,L)) 

RVM1  =  Q(1,L)«C0$(V(NM1,L))-C(2,L)«SIN(V(MM1,1)) 

CVM1  =  Q(  l,U«SIN(V(NMl,l>  )4-G(2»L)«CCS(V(KPltl)  I 
C C  M i  =  ( RVM-RVM 1}*C0S(B*C«5) 

C  5  Ml ) «  (RVM*RVMl)«SIN(B*C.5> 

C<  MN ) *  (CVM-CVM1)«C0SCB*C.5) 

800  C(M2)a  <CVM*CVMl)*SIN(8*C.5) 

C  NOW  THE  MATRIX  IS  FILLED  UP  AND  THE  SOLVING  CF  THE 
C  EQUATIONS  WILL  START 

IF(L-l)  7CC,6CC,7C0 

600  CALL  $CLVE(2«N,A,C|1»C.C1,5«X,IT) 


o  r>  r> 
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GO  TO  52 

700  CALL  SOLVE! 2»N» A,C»  2,0. Cl ,5  » >  »I  T) 

52  N2*  2»N 

00  51  KM«1,N2 
51  ll!KM,L  )*X  ( KM ) 

50  CONT INUE 

NOW  THE  RERAOIATION  PATTERN  WILL  BE  CALCUIATEC 

IF t IT)  S4, 55,55 
55  KLM  *  0 

210  IF ( KLM )  230,220,230 
220  FI»LFI+18C 
KT*  SC 
GO  TO  240 
230  FI»LFI 
KT*1CC 

240  OCALL  PATT(KLM,KT,K, I ,N , J ,0 ,R , W»ZC ,FII ,LTETA,F I ,TETAI tH,U,G»T»F,CLE 
1) 

KLM-KLM+l 

IF(KLM-l)  281,210,281 
281  WRITE  ( 6, 3CC ) 

WRITE  (6,71)  N.PAGINA 
WRITE  (6,72) 

WRITE  (6,74) 

IFIJ.EQ.l)  WRITEI 6, 76) 

IFIBW.EQ.C)  GC  TO  31C 
WRITEI 6, 2CC)  FAKTOR 
IF(MDIP.EQ.O)  WRITEI6.333) 

IF(MO.EO.C)  WRITE<6,444) 

IF(MR.EQ.C)  WRITE(6,555) 

IFIMB.EQ.O)  WR|TE(6,666) 

310  IF(MDB.NE.C)  GO  TO  120 
CO  111  MK  *  1,15 
CO  111  KK  *  1»K 
111  G(MK,KK)  *  OBIG(MK,KK) ) 

120  AA  “  B/6®2£31E531 
DPI*D/6. 28316531 
77  WRITE16, 177)  DPI 
WRITE! 6, 777)  RAO 
WRITEI6, 686)  OLE 
WRITE  <6,175)  AA 
HH=R/6. 28318531 
IF! J-l )  180,178,180 
178  WRITE  (6, 17S )  HH 
180  WRITE! 6, 181)  ZO 

IFIMOB.NE.C)  WRITE !6,3C) 

WRITE  (6,78)  (TETAI(K) ,P*1,K) 

WRITE  (6,7S)  (PHI1 (M) ,M*1»K) 

WRITE  (6,80  <VV!M),K=1,K) 

WRITE  (6,81) 

C 

00  84  MK*1,  IS 

84  WRITE  (6, SC)  T(MK)»F(MK),(GIWK»M) »M*l »K) 

IF(MDB.NE.2)  GO  TO  S6 
WRITE(6,3C) 
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CO  112  MK  -  1,19 
CO  112  KK  >  1,K 

112  6(MK,KKI  •  DB! G!MK,KK) ) 

CO  113  MK  >  1,19 

113  WRITE! 6»9C 1  T ( MK ) »F ! KK ) , (G I MK ,M) , N-l , K ) 

CO  TO  96 

96  WRITE! 6, 95  ) 

96  O-CPI 
B-AA 
R-FH 

IFtBW.EQ. 1 )  GO  TO  504 
IFJCI.NE.e.ANO.I.LT.IMX)  GC  TC  501 
Dl-0 
I-IST 

IF( DC .NE • 0 • .ANO.D.  LT.DMX )  GO  TO  510 
CD-C. 

C-CST 

IFJCC8.NE.0..AND.B.IT.BMX)  GC  TC  512 
COB-O. 

B-BST 

IF!OZ.NE.C..AND.ZO.LT.ZMX)  GC  TO  513 
CZ»0. 

ZO-ZST 

1FIDR.N6.C..AND.R.LT.RMX)  GO  TC  511 
CR-'O  . 

R-RST 

IF!  CD’-E.NE  .0.  .  AND.DLE.  LT.DLKX )  GC  TO  514 
DDLE-O. 

CLE*OLST 

C 

504  IF(BW.EQ.O)  GO  TO  505 
C 

IFIMO.EQ.CI  0»0*FAKTCR 
IF1MR.EQ.C)  R*R*FAKTCR 
IFt MB.EQ.G )  B-B*FAKTCR 
OLE«OL6*FAKT1 
FAKTOR*  FAKTOR+RATIO 
IFIFAKT0R.LT.1.C+0IFF)  GC  TC  502 

505  IF(MA)  93,100,93 
93  CONTINUE 

300  FORMAT  (1H1I 

71  OFORMAT! 31H  SQUARE  VAN  ATTA  REFUECTCR  WITH,I3,10H  ELEMENTS ,50X, 4FP 
1AGE, 14) 

72  OFORMAT! 1C2H-CALCLLATED  RERACI ATICN  PATTERN  *G«  IN  THE  DIRECTION  Cl 
1VEN  8V  THE  ANGLES  TETA  ANO  FI  FCR  VARICUS  CASES) 

74  0F0RMATC67H  OF  INCIOENCE  ANC  POLARIZATION  -GIVEN  BY  THE  ANGLES  TETA 
1I.FII  AND  V) 

76  FORMAT 1 48HCHERE  THE  CONOLCTING  PLATE  IS  TAKEN  INTO  ACCOUNT ) 

200  0F0RKAT(45H0BANCWI0TH  CALCULATIONS  FCR  LAPBOA/LAMBCAO  «,F8.3,4H  AN 
1C) 

333  OFORMAT ( 62H  NO  VARIATION  CF  PHYSICAL  DIPCLE  LENGTH  ! LENGTH*  0.5«  LA 
1MB0AC) ) 

444  FORMAT ( 5CH  NO  VARIATION  CF  PHYSICAL  DISTANCE  BETWEEN  CIPCLES) 

555  FORMAT ( 43H  NO  VARIATION  CF  PHYSICAL  HEIGHT  CVER  PLATS) 

666  FORMAT ( 4 1H  NO  VARIATION  CF  PHYSICAL  LENGTH  r.®  LINES) 

177  FORMAT  (  27H0D I  STANCE  BETWEEN  ELEMENTS- ,F10.C.  „ ,  i  3H  WAVELENGTHS) 


o  r>  r> 
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777  FORMAT ( 27H  OIPOLE  RADIUS  «  ,Fll,3,12H  WAVELENGTHS  l 

888  FORMAT ( 27H  OIPOLE  LENGTH  «  ,Fn.3e12H  WAVELENGTHS) 

175  FORMAT  l  29H  LENGTH  OF  TRAN$MI$S1CNLINES-»F8.2,13H  WAVELENGTHS) 

179  FORMAT ( 33H  DISTANCE  FROM  ELEMENTS  TC  PLATE- ,F4.2,i3H  WAVELENGTHS) 
181  0 FORMA"" ( 46H0CHARACTERI STIC  IMPEDANCE  CF  TRANSMISSION  LINES- , F5 .2, 6H 
1  OHMS) 

30  FORMAT! A8HGRERADI ATION  PATTERN  VALLES  MEASUREC  IN  CECIEELS) 

78  FORMAT! 9H-TETA I*  ,2X,10F9.1) 

79  FORMAT ( 9H  FII  ■  ,2X,1CF9.1) 

80  FORMAT ! 9H  V  «  ,2X,1CF9.1> 

81  FORMAT! 11HCTETA  FI) 

90  FORMAT ( 1H  ,F4. 1,F6.0, 1CF9. 2 1 

95  FORMAT! A2H1»#*THE  SYSTEM  OF  EQUATIONS  IS  SINGULAR*-*) 

STOP 

END 


AIBFTC  DEC  I  OECK 


REAL  FUNCTION  OBtX) 

«  DB  *  CONVERTS  THE  RERADIATEC  ENERGY  INTC  DECIBELS 

DB  *  -9CC.C 
XI  *  1C0CCCCC.*X 
KX  *  IFIXIXI) 

IFIKX.NE.Q)  GO  TO  1 
RETURN 

1  CB  «  10«*ALQG1C( X) 

RETURN 

ENC 


uuu 
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8I8FTC  SU61  CECK 

SUB1  -  EFN  SOURCE  STATEMENT  -  IFNtS)  - 


SUBROUTINE  MAT  1 t LE, l , N , J, B , IB ,A  ,H) 

•MAT1*  COMPUTES  THE  FIRST  PART  CF  THE  MATRIX  ECUATICA 

DIMENSION  IB t ( ) * At  72,72) «H t3£,6) 

MM»( N-LE 1/2 
CO  32  M«1,MM 
CO  30  L*  1 » I 
MX»M-L*I 

IFllCC*MX-«)  27,30, 3C 
|  27  IFt  lCO*tMX-LE)450*I-5)  28,29,29 

\  28  LL*N-2*M42 

|  GO  TO  31 

I  29  lL-2«MX4N~2*l*t+2»! 

*  GO  TO  31 

30  CONTINUE 
C 

31  Y*tl  1,31-FLOATt  J)*Htl,5)-HtU»3J*FtOATI  J)«HCLL,5I»«CCSCe«0.5) 
A(M,2*H-l)»  Y 

MN*M4N 

At  MN»  2*M  )  ■  Y 
NMl«2»tN—M )♦! 

NM2«2«CN-M)42 
MLl-P4tN4LE)/2 
ML2«M4(3*N4LEJ/2 
AtM.NMl)-  -Y 
Ai HN»NM2  )*  -Y 
C 

0Y»-(H(l,4)-FLCAT(J)»Htlf6)>H(Ll,A)4FLCATU)«H(U,6))»CCSie«0.5) 
1  4S  IN ( B*C  ,5 ) 

At M, 2*H ) *  Y 
At MN,NM1  )*  Y 
AtM,NM2l*  -Y 
AtMN,2*M-ll«  -Y 

Y*(H( l, 3 1-FLOAT  tJ}*H(l,5)4Ht LL»3  J-FLCAT tJ)iHtLL»5))*SINIB«0*5) 
At  ML  1*  2»M-1 ) »  Y 
At  ML 1,NM 1 ) ■  Y 
At  ML2, 2*M ) ■  V 
AtML2,NM2J«  Y 

OY— tHtl,4)-FL0ATCJ)«Ht  l,€)4HtLL,4)-FLCATtJ>«HlLL,6)MSlNte«0.5) 
i  -CCS (8*0.51 
AIHL1,2*N )“  Y 
At  ML  I,  NM2 )  *  Y 
At ML2, 2*M-  1 ) »  *Y 
At  ML2»NM 1 ) *  -Y 


C 
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SUB  1  -  EFN  SCLRCE  STATEMENT  -  IFMS)  - 

3 2  CONTINUE 

M»CN-3*LE)/2 
LL  =  ( N-LE 1/2 
CO  4C  Ma  1  *  f*M 
ml*mi 

IF(LL-ML)  SSS,fe?,88£ 

H80  CO  4C  L=Ml,Ll 
K  1  =  C 
L2=  L-M 

CO  35  Ml*M,L,L2 
K 1 -K 1 4 1 
CO  34  L 1  =  1 ,  L  L 
^X=N-  1-L  1*1 
IF(2*^X-1)  3  2  »  ?  4  »  3  4 
i?  IB(  K 1 J  «MX ♦  1- 1 
K1=KH  I 
I8(K1)=L1-1 
GO  TO  35 
:•'=  CONTINUE 
3:  CONTINUE 

IB(  5  1  ■  I-  1-  IB  (  2  ) 

10(6  )  =  I-l-18(4) 

IF(  10(  1)- 18(3))  37,37,36 
36  LL  Is*  I*(  1B(4)-IB(2)  )- 1 8  ( 3 )  ♦  I B  ( 1 )  *1 
GO  TO  113 

U  LL  1*  I*(  JB(4)-IB(2))  +  IB(3)-I8(im 
1)3  IF(  I B  (  1)-IB(5)  )  39,39,36 
36  LL2=l*l IB(6)-1B(2))-IB(5)4IB(1)+1 
GO  TO  114 

39  LL2-  I*(  IB(6)-IBC2))  +  rB(5)-IB(l)4l 

)iH  OY«(H(LL1,3)-FLOAT( J)*H(LL1,5)-H(LL2,3)4FLCAT(J)«H(LL2,5)  ) 

1  »COS ( B*C • 5 ) 

A(L,2»«-1)*Y 
At  M, 2*L- 1 )  =Y 
NL=N+L 
nm=n+m 

NM 1=  2*( N-M ) * l 
NL  1=2*(N-L )  41 
NM2=Niy  l4l 
NL2=NLl4l 
LEl*L+(N4LE)/2 
MEl»M4(N4LE)/2 
LE2*L4(3*N4LE)/2 
ME2=M4(3*N4LE)/2 
A(NM,2*L)«  Y 
A(  NL , 2*M ) *  Y 
A(L,N*U)  =  -  Y 
A(M,NL1)»  -Y 
A(NL,NM2)=  -Y 
A( NK , NL2 ) *  -Y 

OY *- (  H(LL  i,4)-FLOAT( J ) *H ( l LI , 6 )-H ( t 12 , 4 ) *FLC AT ( J)*N(L12,6) ) 
1  *COS ( B*C  *  5 ) 

AIL, 2*M)»Y 
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A(M,2*L)«Y 
AINL.NMU-  Y 
A{ NMfNLl  )•  Y 
A(L,N!*2>»  -Y 
AIM.NL2I-  -Y 
A( NL» 2*M— 1 )»  -Y 
A( NM»  2*L-1) *  -Y 

0Y« ( H( LL 1  *  2 I-FLCAT (J)»H(LL1,5)+H(L12»3) -FLCAT (J)*H(LL2t5) ) 

1  *S  IN( B*C.5 ) 

A(LE1,2*M-1)«  Y 
A(ME1,2*L-1)»  Y 
A(LEl.NMl)*  Y 
A<ME1,NL1>»  Y 
A(LE2»2*M)«  Y 
A<HE2,2*U*  Y 
A(LE2,NM2)»  Y 
A( ME 2»NL  2 ) *  Y 

OY«-(H(LLl.4)-FLOATU)»H|LLl*6HH(U2.4»-FlCAT(J|«h(LL2.6)  ) 
1  *S IN  I B  *C  .5 ) 

A(LEl,2»M)«  Y 
A( ME li 2*L ) »  Y 
A(LE1,NM2>«  Y 
A(NE1,NL2)«  Y 
A ( LE2? 2»M- 1 1 «  -Y 
A(ME2,2»L-1)»  -Y 
A(LE2,NM1)«  -Y 
A(  ME  2»NL 1  )*  -Y 
AO  CONTINUE 
999  CONTINUE 
RETURN 
ENC 
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SIBFTC  SUB2  DECK 

SU82  -  EFN  SOURCE  STATEMENT  -  IFNCSJ 


SUBROUTINE  MA72iLE,I,N,J«B,A,H) 

•  MAT 2*  COMPUTES  THE  SECCNO  PART  CF  THE  MATRIX  ECUATICN 

C IMENS ION  A( 72 » 72 ) ,H I  36 , ft ) 

Y*(M  1 »  3 )— FLOAT!  J)«H(1»5))*CCSIB) 

NE*(N*LE)/2 
NE1-I 3»N*LE)/2 
A ( NE « N  )  *  Y 
A  ( NE 1 , N+ 1 1  *  Y 

A(NE'N*1)*  -IH«l,A)-FLOAT(  J)«H  11,6)  >«CCS  IB) ♦SIMP) 

A(  NE 1*  N  )  *  “A(NE, N*l) 

MM-IN-LEJ/2 
CO  A 7  M*1,MM 
LL*(N-LE)/2 
CO  A3  L*l, LL 
MX-M-L*! 

IF(2*MX-l)  A2,A3,A3 
A2  NX-MXM-1 
NY* I- I 
60  TO  AA 
A3  CONTINUE 
AA  IF(2»NX-I)  A6,A6,A5 
A5  LL3*  1*1  U-l  W2-NY)M3-I  )/2+NX 
GO  TO  115 

Aft  U3»  I*C  ( 1—1 )  /  2-NY  )♦  ( 1*I)/2“NX 

115  Y*(H(LL3,3 )- FLOAT  fJ)*H(LL3«5))*CCS(B) 

NM*N+M 
N1MN+1I/2 
N2*( 3*N+ 1 ) tZ 
NM 1*2* ( N-M  )♦  1 
NM2«  NM1+1 
M2*  (3*N41J/24M 
Ml*  (N  +  D/2+M 
AIN1,2*M-1 )*  Y 
A(N1,NM1)«  Y 
A<N2,NM2)«  Y 
A(N2,2»M)«  Y 
C 

Y*-( H(LL3, A ) -FLOAT (J)*H(LL3,6) )*CCS(Bt 

A(N1,2*M)»  Y 

A(NlfNM2)»  Y 

A(N2,2»M-D*  -Y 

AIN2.NM2)*  -Y 
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SUB2  -  EFN  SCARCE  STATEMENT  -  IFMSJ 


Y*2.»(F(LL2*2 )-FLOAT (J)»H(LL3,5))*SIN(B*0.5) 

At  M2  *  N* 1 ) ■  Y 
A(M1,N)«  Y 

At M 1  * N+l ) »  -2.»tH<LL3,4)-FLCATt  J»«HtU3*6))*SIMe»0.5) 
A(M2fN)«  -A ( M 1 »N*1 ) 

A|M,N)>O.C 
At  M*  N* 1 ) *C  »C 
A(NMtN)>  C.C 
A7  At  NM , N+l ) ■  O.C 
RETURN 
ENC 


o  n  n 


-  32  - 


014C02 


SIBFTC  SUB 3  CECK 

SUB3  -  EFN  SOURCE  STATEMENT  -  IFMS) 


OSUBROUT INE  PAT 7 (KIM, KT , K , I *N , J ,0 , R , W, ZC , F I ! , LTET A , F I *TETA I »H»U, 

1  G, T,F»DLE ) 

*PATT»  COMPUTES  THE  RERACIATICN  PATTERN  CF  THE  REFLECTCR 

DIMENSION  TETAK  1C),H(36  ,6)  ,1(72,10)  .6119,10)  ,T(19)  ,F(  19) 

PO(X)=  SIN(0.5*X»XI*D)/(C.5»X»XI*C) 

MKM-C 

CO  91  NTETA*1C,KT,LTETA 

HKHsMKH+1 

MK=  KLM«  9*MKK 

NX  *  {MXM-ll)«(KLN-l)+  HKP  *  KLH 
IF(KLM  )  260, 25  0,260 
250  MTETAMCC-NTETA 
GO  TO  270 

260  MTETA»NTETA  -  1C 
270  00  69  1 , K 

GR=0. 

G I  *0. 

Y=3,1A159265/160. 

TETA.MTETA 
TETA  ■  TETAI (NX  ) 

T(HK)*TETA 

F(MK)»FI 

0XKl*120.*3.lAl552fc5*(CGS(3.iAl59265*DLE«SIN(TETA*Y)*SIN(FI«V) ) 

1  -COS(  3.14155265*OLE>)  /  (ZC*  SORT  C 1.- CSI MTETA«Y)  *S I M  FI  *Y)  )  »•£ 

2  SIN(3.1A155265*0LE)) 


84 

85 


N2»2«N 

CO  64  L  =  1,N 2,2 

Ll*(L+l)/2  ,  .  . 

C1*H(L1,1)*SIN(  TETA*  Y)*CCS  (FI*Y)  ♦«  Hi,  2)  *SIMTET  l*Y)«SIMFI*Y) 
AA*  U(L,M) 


NL“L ♦ 1 
BB=  UCNL.M) 

GRaGR*AA*COS( C 1 ) ♦BB*SIN  <C 1 ) 

GI*GI-AA*SIN(CI)+BB»C0S(C1) 

CONTINUE 

IF(J-l)  87,85,87 
CO  86  L* 1 , N2, 2 
Ll*(L+l)/2 

OC l*(-( L  1,  1  )*SIM  TETA*Y)*CCS(FI*V)+H(ll,2)*SIN(TETA«Y)  «S!N(FI«Y) 
l  -2.*R*C0S(TETA*Y) 

AA*  UU,M) 


NL  »L ♦ 1 


BB*  U(NL,M) 
AN*N 


)• 
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SUB 3  -  GFN  SCLRCE  STATEMENT  -  IFMS)  - 


CR=GR-AA«CUS»Cn-»B«$I\(Cn 
GI  =  GI+AA»S!N(C1  J-BB«CCS<Cn 
8t  CONTINUE 
X  I  «=  I 

AR*COS(R*(COS(TCTAltP)*Y)+CCJ(TFTA*Y>)) 
AI=SIN(R«(CC£(TCTA!ty)*Y)+CC'(TETA«Y) ) » 

Al.  =  SlN(TeTAI(M*Y>*CCS<EIl«Y>  +  SIMTeiA*Y)«CCS(FI*Y) 

BF*SIN{TETAI(^J*Y)»SIN(FII*Yl*SIN(TETA»Y)*SIMFI*Yl 

AL1=AL*1CCCCC. 

flEl=BE«lCCCCC. 

IF (  IFIX(ALl)  )  29C,291,29C 

291  SAL5'  1  • 

GO  TO  292 
29C  SAL  =  PO(AL) 

292  IF C  IFIXCBEI) )  294,292,294 

293  SBE= 1« 

GO  TO  29 £ 

294  SBE  =  POtBE) 

295  XK  2=  C,5«SAL*SBE*AN*C**2 

TP=COS(W«Y)*CCS<FII*Y)*SIMfc»Y>*CCS{TGTAUw>«Y)*$IMFI  I«Y  I 
TT*-CO$(k*Y)*$IN(FII*Y)  +  <lM  k*Y)«CCS  (TETAI  (F)«YJ»CCS(FII«Y) 
FT=-TP*CC$(FI»Y)*CC$(  TETA«Y)*TT«SIMFI*Y)*CCS(T£TA*Y) 
FP*TP*SIN(FI*Y)+7T*CCS(FI»Y) 

OG(PK,*»*  (  (XK1#GR«S!N(FI*Y»  •CCS(TETA«Y)  ♦  XK2«AR«FTM*2 

1  +<XK1*GI*SIN(F1*Y>*CCS(TETA*Y)  ♦  XK? • A  I *FT > ««2 

2  +(XKl*GR*CCS(FJ«Y)  ♦  XK2»AR«FP1 »  •  2 

3  ♦(XK1*GI*CCS(F!«Y)  ♦  XK2M I  *FP )  **2  )  /  l  3. 14 159265«« 3  ) 
fi7  IFU-1)  6€,€9,E6 

88  CG(MK,y)  =  XKl»*2*(GR‘*2+Gl»*2)*nSlMFI*Y))««2«(CCS(TETA»Yl)««2 
1  +(CCS( FI*Y  I  )••?)/( 3, 14159265**3  I 

89  CONTINUE 
91  CONTINUE 

RETURN 

ENC 
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tieFTC  IMP Z  CECK 

impz  -  efn  source  statement  -  ifmsi 


SUBROUTINE  BETA ( II , J ,  ZC ,R ,H , FAKTCR ,RAD ,OLE ) 

•  BETA*  BULDS  LP  AN  ARRAY  CF  SELF-  AND  MUTUAL  1 MFECANCES 

CIMENSION  M26,6) 

CALL  SELFICLE/FAKTOR,RAO,RS,XS) 

CO  2C  M* 1* II 
IFIM-1I  17,16,17 

16  H  1,  2  J  ®  RS/ZC 
Ml, A)  =-X</ZC 
CO  TO  18 

17  0Z»  AMPE(CLE/FAKTOR,OLE/FAKTCR, 

1  A8SIC.15915A5*(H(l,l)-H(M,im  , 

1  AB S I C. 1591 549* <H 1 1,2 1-H I M,2) > ) ,0.0,0.0,1,0.001,10) 

MM,  3)«Z/Z0 

0 Z*  AMPE(CLE/FAKT0R,0LE/FAK7CR, 

1  ABS(C.1591549*lHtl»l )— H ( M, 1 )  )  >  , 

1  A0S(C.1551549*(H<1,2)-H(M»2)))*C.O»O«G»O»O.OO1»1O) 

MM.4MZ/20 

19  IF ( J )  19,20,19 

19  0Z=  ampe«cle/faktgr,dle/faktcr, 

1  0. 1591549* SORT! <H I 1, l )-H (M, 1 )> **2+4. ■( R**2 ) > 

1  AB SIC. 1591549* (HI  1, 21-HI M, 2 ))> ,0.0,0.0,1,0.001,10) 

MM,  5  ) *Z/ZC 

0Z=  AMPEI cle/faktor.ole/faktcr, 

1  C. 1591 549* SQRTI IH 1 1 , 1 )-H (M , 1 ) > **2*4. *  I R«*2 >  ) 

l  ABSIC.1591549MHI  1,2) -HI M,2) ) ) ,0.0,0.0,0,0.001,10) 

MM,6)*Z/Z0 

20  CONTINUE 
RETURN 
ENC 


nnn 
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♦ I6FTC  XCAL  CECK 

XCAL  -  EFN  SOURCE  STATEMENT  -  IFNIS)  - 

SUBROUTINE  SELF < H, A, R , X ) 

•SELF*  COMPUTES  THE  SELF  IMPEDANCE  CF  THE  CIFCLES 


COMPLEX  CS1NT *CINSI 
PI*2 •  1A159265 
AA*SQRTI 2. )*2.*PI*A 
FH*P I*F 

CSINT*CIN$I(2«*HH) 

CI2*  C.57721S££5+AL0G(2.*HH)-REAL ICSINT) 

S2*  REAL ICSINT) 

S 12*  A  IMAG ICSINT ) 

CSINT*  CINSI I  A.*l-H) 

CIA*  0. 57?215£65*AL0G ( A. *HH) -REAL ICSINT* 

SA*  REAL ICSINT) 

S  I  A*  AIMAGICSINTJ 

X*  -29. 997925* I S INI2. *HH )» I-C. 577215665+ ALCGIHH/ I AA**2») *2. *C 12 
C  -CIA)  -  COS  1 2. *HH) • 1 2.*SI2-SIA )-2.*SI2) 

R*  29.997525*II2.+2.*CCS:2.*HH))*S2-CCS12.*HH)*SA 
C  -2.*SINI2.*HH)*SI2+SIN12.*HH)*SIA) 

X=X/l SINIHH)**2) 

R*R/ISIN(HH)**2) 

12  CONTINUE 

13  CONTINUE 
RETURN 
END 
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$  IBFTC  CINSI  DECK 

CINSI  -  EFN  SOURCE  STATEMENT  -  XFNIS) 


COMPLEX  FUNCTION  CINSI (X) 

C  CINSI  COMPUTES  AS  ITS  REAL  PART  THE  MODIFIED  CCS IKE  INTEGRAL  AKC  AS 
C  ITS  IMAGINARY  PART  THE  SINE  INTEGRAL. 

IF  (X.GT.1.1)  GO  TO  20 
C  =  1.C 
S=1.0 
TC»1.C 
TS  = 1.0 
Y  =  X«X 
1=0 

10  1  =  1*1 

TC=-Y«TC/FLOAT I 2»I* ( 2*1- 1 ) ) 

TS*-Y*TS /FLOAT 12*I*(2*I+1) ) 

TERMC=TC /FLOAT  C  2* I) 

TERMS*TS /FLOAT  12*1*1) 

C=C*TERMC 
S=S+TERMS 
EC*ABS( TERHC/C ) 

ERRQR=l.CE-8 

IF  I  EC. GT. ERROR)  GO  TO  1C 

CIN*1.0-C 

SI»X*S 

C INS I*CMPLX(C IN t  SI ) 

RETURN 
20  Y=X*X 

F= l Y*I Y*  t  Y*( Y+38.027264 ) *265. 187033 ) *335.677320 ) +38. 102A95 ) / 

C  « X*C  Y»(  Y»(Y»(Y*4C. 021433) *322. 624911 »*57C. 236280)  *157.105422) ) 
G*tY*CY*CY*(Y+42. 2428551*302. 757865)  *352. 018498)*21. 8218991/ 

C  (Y«{Y*{Ys{Y*{Y*48.19c527  } *482. 465984 )*ii 14.978685  1*449 .6903  26  i ) 
S*SIN( X) 

C»CQS<X) 

CIN=0.577215665*ALOG(X)-F*S+G*C 
S  1=1 . 57079633- F#C-G*  S 
CINSI»CMPLX(CIN,SI) 

RETURN 

END 


nonoonn 
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tlEFTC  ZCAL  CECK 

ZCAL  -  EFN  SOURCE  STATEMENT  -  IFMS)  - 


REAL  FUNCTION  AMPE I H 1 ,H 2 , YC , Z C  ,  TF T A , F I , NC , CE L T A , FCRC ) 

■AYPE»  COMPUTES  THE  "LTLAl  lYPECANCCS  GETfeEEN  THE  CIFCLES  EY  AN 
INTECRAT ICN  USING  RCFBERG'S  PETHCC 

THE  REAL  PART  CF  THE  *1 TLAL  I  &PECANCE  FC«  ACM  AN C 
TFC  IMAGINARY  PART  FCR  NC=C 

CIMENSION  TRAPU1) 

STEP*  F2 
X  =  -H2/2* 

Cl*RX(Fl,F2,YC,ZC,TCTA,Fl ,\C,X) 

X=FZ/2 . 

FC *RX { H 1, F 2, YC, ZC.TETA.FI ,NC,X) 

TRAP( l)=tCl*FC )*STEP/2. 

C1=C. 

CO  2C3  K=1,MCRC 
SUM=  C • 

ERRCR=C. 

STEPsSTEP/2. 

M*2#*K 

PV=P-I 

CO  2 C4  L»1.mp,2 
LL=M-L 

XL=FlOAT(LL) 

XM*FLO AT { P ) 

X=XL/XK 

X  =  XM-F2/2.M(l.-X)*F2/2. 

FO=RX < HI, F2, YC.ZC.TETA.FI ,NC,X) 

C2sSUmfFC 

in  ABS(  FC  )-AR<(  5LMn  2C6.2C6.2C5 
205  ERP.CR*ERRCR*$LP-(C2-FC) 

CO  TO  207 

2CC  ERRCR«ERROR+FC- I C2-SUW) 

2C7  SU^=C2 
2CA  CONTINUE 

TRAP (K ♦  ! )«TRAP (K1/2.+ JC2+ERRCR) -STEP 
P=  1 . 

KK=IC-K 

CC  2CP.  LL  *KK ,9,1 
L=-(LL-1C> 

P*P»A. 

2 08  TRAP( L )=(TRAP(L+1>»P-TRAP<L)I/(F-1.) 

C2*  TRAP (  1  > 

IF(K-1>  2C9. 2C9.211 

211  ir(ADS(C2-Cl)-CELTA*ABS(C2))  21C.210.2C9 

209  C 1 =C 2 


3  C  - 


C14CC2 

7CAL 


EFN  SCIRCE  STATEMENT  -  imS) 


203  CONTINUE 


MORC=K-l 

210  0  =  2,>, 97925 /SIN (3.141 
IF(NO-l)  2C 1 , 2CC»  2C1 
2CC  AMPf-  =  -Q*C2 


GOTC  202 

201  AMPE*  0*C2 

202  CONTINUE 
RETURN 
EMC 


59265*H1)/$IIS 


(3. 14159265-H2) 


r»  ^  r> 
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$ I8FTC  INTG  DECK 

INTG  -  EFN  SOURCE  STATEMENT  -  IFMS) 


REAL  FUNCTION  RX(H1 ,H2 , VC , 20 , TETA ,F I ,NC , X) 

*  RX  *  COMPUTES  THE  FUNCTION  TO  BE  INTEGRATED  IN  ROUTINE  AMPE 

b 

B2»F2/2. 

A=CO$<TETA) 

3»SINI TETA} 

C*B*SIN( F ! ) 

C«B»COS(FIi 

E*2.*CQSI  2«14159265»Hll 

F«Z04B1 

G-ZO-61 

$X»X*C 

SY*X»C 

SZ*X#A 

R02=SX«»2  4(YC4SY)**2 
H*Z04SZ 

R»SQRT(R024H**2  J 

FH-F4SZ 

P-G4SZ 

R1«SQRT(RC24HH**2  ) 

R2*SQRTtR024P**2  ) 

IF(NO-I)  214,211,214 

211  AK=SINC6.2831€531»R11/R1 
AL*S INI 6.28316531*R2)/R2 
AM *SIN(6. 2831653 1*R)/R 
GO  TO  212 

214  AK=COS (6. 28318531 *R1)/R1 
AL»COS(6.26316531*R2)/R2 
AM*CQ$( 6. 2831653 l*R)/R 

21 2  CONTINUE 

I F t YO  I  218,215,218 

215  IF C  SX I  218,216,218 

216  IF ( $Y 1  218,217,218 

217  RX«lE*AM-AK-AL)*A*SIN<6.2e3ie53l»(B2-ABSIxn I 
GO  TO  219 

2180RX-I ( AK*hH4Al*P-E*AM*H)*{X*ID*»2}4YO«C4X« (C««2))/RC2 
1  4( E*AM-AK-AL )*A )• SIN (6* 28318531* (B2-ABSIX) )  ) 

219  CONTINUE 
RETURN 
END 
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ilBFTC  SOLV  DECK 

SOLV  -  EFN  SOURCE  STATEMENT  -  IFNCS) 


SUBROUTINE  SOLVE  INN, A,B ,  I N ,E P S,  I TMX  ,X,  IT) 

C 

CSOLVE  LINEAR  EQUATION  SOLVER  WITH  ITERATIVE  IMPROVEMENT  VEPS1CN  IV 
SOLVES  AX*B  WHERE  A  IS  NXN  MATRIX  AND  B  IS  NX1  VECTOR 
IN* 

1  FOR  FIRST  ENTRY 

2  FOR  SUBSEQUENT  ENTRIES  WITH  NEW  B 

3  TO  RESTORE  A  AND  B 

EPS  AND  UMAX  ARE  PARAMETERS  IN  THE  ITERATION 
IT* 

-1  IF  A  IS  SINGULAR 
0  IF  NOT  CONVERGENT 
NUMBER  OF  ITERATIONS  IF  CONVERGENT 
CALLS  MAP  SUBROUTINES  ILCG2 ,OCT, SOOT  AND  CAD 

TO  MODIFY  DIMENSIONS*  CHANGE  THE  NEXT  3  (NCT  2  BUT  3)  CAROS. 

OCIMENSION  A (  72, 72) ,  B  172)  ,  X  ( 721  ,DX  172)  >R  (72  )  »Z  ( 72)  »  RMI72  J  .  IRP 1 72  ) . 
1AA( 72*72) 

MA*72 

MA  MUST  *  DECLARED  OIMENSICN  CF  SYSTEM 
EQUIVALENCE (R*OX ) 

GO  TO  ( ICCC,2CCC,3CCC) , IN 
1000  N=NN 

NM1-N-1 
NPl*N-»  1 

EQUILIBRATION 

CO  S 1C  1*1, N 

KT0P»IL0G2(AtI*in 
DO  5C3  J-2,N 

503  KTOP*MAXC(KTOP, I  LOG 2 (A ( I  *  J) I ) 

RM(  n»2.C**l-KTOP) 

DO  5C<5  J«1,N 

509  A(I,J)*AII,J)*RM(I) 

510  CONTINUE 

SAVE  EQUILIBRATED  DATA 

DO  54C  1*1, N 
00  5 AC  J*1,N 
548  AA( I  * J )*A( I »  J ) 

GAUSSIAN  ELIMINATION  WITH  PARTIAL  PIVCTING 
CO  59  M*l,NMl 


nnn  n  o  r>  r» 
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SOL V  -  EFN  SCLRCE  STATEMENT  -  IFN(S)  - 

TCP-ABS  <A(P,P)) 

IMAX-M 
CO  12  i*P,N 

I F  {  TCP-ABS  1AU  ,F)  1110,12,12 
10  70P*ABS  ( A( I ,K) ) 

I MAX* I 

12  CONTINUE 

IF  (  TOP ) 14  *  13  < 14 

13  IT— 1 

C  *$ INCULAR* 

RETURN 

14  IRP ( M ) • JNAX 

23  IF(  IMAX-PI29»29,24 

24  CC  25  J*1,N 

T£MP*A(H, J) 

A  CM, J ) *AC I MAX  *  J) 

25  ACIMAX,J)«TEPP 

29  MP1-H+1 

CC  33  I*PP1,N 

EM*ACI,M)/A<M,M> 

AC  1 ,M)*EM 
I F ( EM ) 3 1, 33, 3 1 

31  CC  32  J*PPi,N 

32  All,J)-AlI,J)-AtK,J)»EM 

33  CONTINUE 
99  CONTINUE 

IRP(N)*N 

IF  (AIN, NH12C, 113,120 
113  IT  — 1 
RETURN 

120  CONTINUE 

STORAGE  FOR  A  NON  CONTAINS  TRIANGULAR  L  AND  U  SC  TEAT  CL*I)«U»A 

CUPLICATE  INTERCHANGES  IN  OATA 

CO  229  1*1, N 
IP»IRP(  I J 

IF ( I- IP ) 221 , 229, 22 1 

221  CO  222  J*1,N 
TEMP«AAU,JJ 
AA( I , J)*AA( IP, JJ 

222  AA ( IP, J I *TEHP 
229  CONTINUE 

PROCESS  RIGHT  HANO  SIOE 

2000  CONTINUE 

CO  6C1  1*1, N 
601  B(I)*8m*RPCn 

CO  609  1*1, NP1 
IP* IRP (  I  ) 

TEMP*B(I  ) 

si  n*Bi  ip) 

B( IP )=TEHP 
609  CONTINUE 
C 
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SOL V  -  EFN  SOURCE  STATEMENT  - 

t 

C  SOLVE  FOR  FIRST  APPROXIMATION  TC  X 

C 

199  00  2C0  1*1 .N 

2C0  ZU)=-SD07(  I-1,A(  1 ,1)  ,MA,2m  ,1,-8111) 

CO  201  K* 1 ,N 
I*NP1-K 

201  xm*-SOOTtN-IfA(I.I«l)tMA'X(I*l)  ,1.-Z1I) 1/AII *1) 

C 

C  ITERATIVE  IMPROVEMENT 

C 

IF<  UMAX )37C,370,3C0 
300  TOP-O.O 

DO  302  1*1, N 

303  T0P«AMAX1(T0P,A8S(X(I>)> 

EPSX*EPS«TOP 
DO  269  IT  l.ITMAX 
C  FIND  RESIDUALS 

DO  319  1*1, N 

319  RJ  I)=-D0T(N,AA(I»1),KA»X(1)  ,1,-BU)) 

C  FIND  INCREMENT 

00  329  1*1, N 

329  2(I)t-S00T(M,A(l,ll,FA,mM,-RU)) 

00  339  K*1,N 
I*NP1-K 

339  0X(  I  I*-SD0T(N-I,AU,I  +  1)  ,MA,DX(I*1)  ,1,-ZIIII/A!  1,11 

C  INCREMENT  AND  TEST  CONVERGENCE 

TOP-C.C 
00  392  1*1, N 
TEMP-XI II 

X<  I)*CAO(Xm,DX(n) 

CELX=ABS  ( X ( X  * -TEMP) 

TOP»AMAX1(TCP,OELX) 

392  CONTINUE 

IF (TOP-EPSX 1381, 38 1,369 

369  CONTINUE 

370  IT*0 
381  RETURN 

RESTORE  A  ANO  8 

3000  CONTINUE 

DO  7C9  K*1,N 
I*NP1-K 
IP»IRPI I | 

IF ( I-1PI7CI, 709,701 

701  TEMP*BU) 

B( 1  )*B( IP) 

6< IPI-TEMP 
00  702  J*1,N 

TEMP*AA ( I , J I 
AA(  I,J)-AA( IP, J) 

702  AA( IP, J ) “TEMP 
709  CONTINUE 

CO  729  1*1, N 

Bm-em/RMin 
00  729  J*1,N 


o  <_>  o 
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SOLV  -  EFN  SCLRCE  STATEMENT 

A(  I,  J)«AA{  I,  Jl/RNII  ) 

729  CONTINUE 
RETURN 
ENC 


*  $  I BFTC  LOG  DECK 


INTEGER  FUNCTION  ILOG2U) 

■ ILOG2*  ROUTINE  POR  USE  H  TH  ROUTINE  SCLVE 
ILC62*C 

IF  (Z.NE.C.)  GC  TO  1 
RETURN 

1  IL0G2*A(NTO.2  2  2C»ALCGlC(ABS(Zm 
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Fig t. Coordinate  system  for  arbitrary  Van  Atta  array* 


Fig 2.  The  Van  Atta  principle  (retrodirective  effect). 
The  paths  A  BCD  and  A'B’C’D'  are  equal. 


Fig 3.  The  specular  reflection  (mirror  effect). 
The  paths  ABC  and  A'B'C*  are  equal. 
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Fig. 5.  Normalized  bock- scattering  cross  section 
of  16 element  square  Van  Atta  reflector 
with  conducting  plate. 
asO.AVX  ,  d=0.6Xe  hsO,25X .  Z#=73ohms, 
*s(f.  v=90! 


rin.tf.  itarmolteed  back-ucrtterirtg  cross  section 
as  a  function  of  the  length  of  the  dipoles. 
N*?$  elements ,  a  s  0.41X ,  d  «  0.6X ,  h  3  0. 25  X 
Z= 73  ohms  .dipole  radius *0.01 5X 
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Fig.7. 16  element  square  Van  Alta  reflector  without  plate, 
a  =0.41  X.  ,  d=0,6X*  ,  Z.  =  73ohms 
dipole  length  =  0,5X#  .dipole  radius  =  0.015A, 
fs  =  0°,  v  =  90° 
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a.  Bandwith  properties. 


b. Measure  M  of  the  retrodirective  effect. 


FigA.16  element  square  Van  Atta  reflector  with  plate, 
a  =0,41  A.  ,  d=0,6X.  ,  h  =  0,2  5\.  ,  Z.  =73 ohms 
dipole  length  =  0,5A#  .dipole  radius  =  0,01 5 A, 
fi  =  0°,  v  =  90° 
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